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Abstract 



The known logarithmic extensions of conformal and Schrodinger-invariance assume 
translation-invariance in their spatial and temporal coordinates. Therefore, they can- 
not be applied directly to slow far-from-equilibrium relaxations, where time-translation- 
invariance no longer holds. Here, the logarithmic extension of ageing-invariance, that is lo- 
cal dynamical scaling without the assumption of time-translation-invariance, is presented. 
Co-variant two-point functions are derived. Their form is compared to transfer-matrix 
renormalisation group data for the two-time autoresponse function of the ID critical con- 
tact process, which is in the directed percolation universality class. 
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1 Motivation and background 



Dynamical scaling naturally arises in various many-body systems far from equilibrium. A 
paradigmatic example are ageing phenomena, which may arise in systems quenched, from some 
initial state, either (i) into a coexistence phase with more than one stable equilibrium state or 
else (ii) onto a critical point of the stationary state [TOl [151 139] . From a phenomenological point 
of view, ageing can be defined through the properties of (i) slow, non-exponential relaxation, 
(ii) breaking of time-translation-invariance and (iii) dynamical scaling. Drawing on the analogy 
with equilibrium critical phenomena, where scale-invariance can under rather weak conditions 
be extended to conformal invariance [661 E] > m recent years it has been attempted to carry out 
an analogous extension of simple dynamical scaling, characterised by a dynamical exponent z, 
to a new form of local scale-invariance (lsi). One of the most simple predictions of that theory 
is the form of the linear two-time autoresponse function [33, 35, 36J 
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which measures the linear response of the order-parameter 0(i, r) with respect to its canonically 
conjugated external field h(s, r). The autoresponse exponent \r and the ageing exponents a, a' 
are universal non-equilibrium exponents^] The causality condition t > s is explicitly included. 
The foundations and extensive tests of f 1 1.11) are reviewed in detail in 



In the case of a degenerate vacuum state, conformal invariance (of equilibrium phase tran- 
sitions) can be generalised to logarithmic conformal invariance [221 EZ], with interesting 
applications to disordered systems [13], percolation [211 EE] or sand-pile models [65]. For re- 
views, see [20, 24|. Here, we shall be interested in possible logarithmic extensions of local 
scale-invariance and in the corresponding generalisations of (II. ip . 

Logarithmic conformal invariance in 2D can be heuristically introduced [261 EZ] by replacing, 
in the left-handed chiral conformal generators £ n = —z n+l d z — (n+ l)z n A, the conformal weight 
A by a matrix. Non-trivial results are only obtained if that matrix has a Jordan form, so that 
one writes, in the most simple case 

£ n = -z n+1 d z -(n + l)z n ( * } ) (1.2) 



x A , 

Then the quasi-primary scaling operators on which the £ n act have two components, which we 
shall denote as $ := ( , ] . The generators (11.21) satisfy the commutation relations [£ n ,£ m ] = 



{n — m)£ n+m with n, m 6 Z. Similarly, the right-handed generators £ n are obtained by replacing 
z i—)- z and A i— > A. A simple example of an invariant equation can be written as = 0, with 
the Schrodinger operator 

S:=(° 9 t) (1.8) 

and because of [«S,^ n ] = — (n + l)z n S — {n + l)nz n+l ^ J d s , one has a dynamic symmetry 
of S^fr = 0, if the conformal weights A = A = are chosen. 



1 In magnets, mean- field theory suggests that generically a — a' for quenches to T < T c and a ^ a' for T = T c 
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Of particular importance are the consequences for the form of the two-point functions of 
quasi-primary operators, for which only co-variance under the finite-dimensional sub-algebra 
(4fci,o) — s[(2,M.) is needed [261 EZ] (we suppress the dependence on z±, but see |T6] ) - Set 

F '■— (0i(^i)02(^2)} , G := (Mzi)Mz2)) , H := (^1(^2(^2)) (1.4) 

Translation-invariance implies that F = F(z),G = G(z) and H = H(z) with z = Z\ — z<i- 
Combination of dilation- and special co- variance applied to F,G leads to A := Ai = A2 and 
F(z) = 0. Finally, consideration of H(z) leads to 

G{z) =G(-z) = G \z\~ 2A , H{z) = (H -2G ln\z\)\z\- 2A (1.5) 

where Go, Ho are normalisation constants. We emphasise here the symmetric form of the two- 
point functions, which does follow from the three co-variance conditions. 

Recently, 'non-relativistic' versions of logarithmic conformal invariance have been studied 
[4Tj . Besides the consideration of dynamics in statistical physics referred to above, such studies 
can also be motivated from the analysis of dynamical symmetries in non-linear hydrodynamical 
equations [6U EH H21 [291 EQ], or from studies of non-relativistic versions of the AdS/CFT 
correspondence [52], El [69l EH [22l HTJ [28] . Two distinct non-semi-simple Lie algebras have been 
considered: 



1. the Schrddinger algebra stf)(d), identified in 1881 by Lie as maximal dynamical symmetry 
of the free diffusion equation in d = 1 dimensions. Jacobi had observed already in the 
1840s that the elements of scl)(d) generate dynamical symmetries of free motion. We 
write the generators compactly as follows 

y«) = -t m+1 / 2 dj -(m+ \)t m -V 2 rj 

M n = -t n M (1.6) 
= -t n ( rj d k -r k d,) 

where Ai is a dimensionful constant, x a scaling dimension, dj = d/drj and j,k = 
l,...,d. Then sct)(d) = (X ±lfi , Y^ /2 , M , B$'%, k =1 d is a dynamical symmetry of the 
free Schrddinger equation S(f) = (2Md t -V 2 r )(f) = 0, provided x = d/2, see [i5l [271 [571 113] . 
and also of Euler's hydrodynamical equations [SI]. An infinite- dimensional extension is 

(X n ,Y m , M n , R n )neZ,meZ+|,j,fc=l,...,d [S]- 

2. The Schrddinger algebra is not the non-relativistic limit of the conformal algebra. Rather, 
from the corresponding contraction one obtains the conformal Galilei algebra CGA(d) [30] . 
which was re-discovered independently several times afterwards [321 [561 ESI El EQ]. The 
generators may be written as follows [H] 

X n = -t n+1 d t - (n + l)t n r ■ V r - n(n + l)r _1 7 • r - ar(n + l)t n 

= - (« + i)t» 7i (1.7) 
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where 7 = (71, ...,7d) is a vector of dimensionful constants and x is again a scaling 
dimension. The algebra CGA(d) = (X±i t o, Y± x , R^^) j,k=i,...,d does (conditional) 
dynamical symmetry in certain non-linear systems, distinct from the equations of non- 
relativistic incompressible fluid dynamics [73j [H]o The infinite-dimensional extension 

\ A Ti) 1 n , -Tin /? 



is straightforward. 



For both algebras sct)(d) and CGA(d), the non- vanishing commutators are given by 

[X n ,X n ,\ = (n-n')X n+n , , [X n ,Y®] = Q -™) Y^ m , >yW] = ^yW-jMyW (1.8) 

where the dynamical exponent z = 2 for the representation (11.61) and z = 1 for the representa- 
tion 1 11. 7ft . For the Schrodinger algebra, one has in addition \Y$,y!$ /2 ] = 5 j ' k M . 

The algebras sct)(d) and CGA(cf) arise, besides the conformal algebra, as the only possi- 
ble finite-dimensional Lie algebras in two classification schemes of non-relativistic space-time 
transformations, with a fixed dynamical exponent z, namely: (i) either as generalised conformal 
transformations jTTj or (ii) as local scale-transformations which are conformal in time [34J. 

Now, using the same heuristic device as for logarithmic conformal invariance and replacing 
in the generators X n in (jl.6|1.7p the scaling dimension by a Jordan matrix 



x 1 — y 



x 1 
x 



1.9) 



both logarithmic Schrddinger-invariance and logarithmic conformal galilean invariance can be 
defined [4Tj . Adapting the definition (II. 4p . one now has F = F(t,r), G = G(t,r) and H = 
H(t, r), with t := t\—tz and r := ri~r2 because of temporal and spatial translation-invariance. 
Since the conformal properties involve the time coordinate only, the practical calculation is 
analogous to the one of logarithmic conformal invariance outlined above (alternatively, one 
may use the formalism of nilpotent variables [HU E]). In particular, one obtains x := x\ = 22 
and F — 0. Generalising the results of Hosseiny and Rouhani j4T] to d spatial dimensions, 
the non-vanishing two-point functions read as follows: for the case of logarithmic Schrodinger 
invariance 

M r 2 ' 

(tfo-GoMWexP " 



G = G \t\ x exp 



subject to the constraint [6J M. : 
invariance 



Mr 2 

YT 



H 



2 t 



1.10) 



-M 2 E For the case of logarithmic conformal galilean 



G — Gq \t\ 



cxp 



t 



H =(H - 2G In \t\) \t\~ 2x exp 



, 7 ■ r 

t 



1.11) 



2 The generator Xq leads to the space-time dilatations t \-> X z t, r Ar, where the dynamical exponent z takes 
the value z = 2 for the representation (jl.6j) of sc\j(d) and z = 1 for the representation (|1.7j) of CGA(rf). We point 
out that there exist representations of CGA(ei) with z = 2 [3S]. From this, one can show that age(l) C CGA(l) 
as well. 

3 In order to keep the physical convention of non-negative masses M. > 0, one may introduce a 'complex 
conjugate' <f>* to the scaling field 4>, with M* = —M.. In dynamics, co- variant two-point functions are interpreted 

as response functions, written as R(t, s) — (<p(t)cfr(s) ) in the context of Janssen-de Dominicis theory, where the 



response field <j> has a mass Ai = —Ai, see e.g. [15J|39] for details. 

Furthermore, the physical relevant equations are stochastic Langevin equations, whose noise terms do break any 
interesting extended dynamical scale-invariance. However, one may identify a 'deterministic part' which may 
be Schrodinger- invariant, such that the predictions (jl.lOp remain valid even in the presence of noise |64) . This 
was rediscovered recently under name of 'time-dependent deformation of Schrodinger geometry' [5 5) . 
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together with the constraint 7 := r y 1 = 7 2 . Here, Go, H are again normalisation constantso 

From the comparison of the results (11.10111.111) with the form (1 1.6ft of logarithmic conformal 
invariance, we see that logarithmic corrections to scaling are systematically present. As we shall 
show, this feature is a consequence of the assumption of time-translation-invariance, since the 
time-translation operator X_i = —d t is contained in both algebras. On the other hand, from 
the point of view of non-equilibrium statistical physics, neither the Schrodinger nor the confor- 
mal Galilei algebra is a satisfactory choice for a dynamical symmetry, since time-translation- 
invariance can only hold true at a stationary state and hence eqs. f ll.6lll.7p can only be valid in 
situations such as equilibrium critical dynamics. For non-equilibrium systems, it is more natural 
to leave out time-translations from the algebra altogether. An enormous variety of physical sit- 
uations with a natural dynamical scaling is known to exist, although the associated stationary 
state(s), towards which the system is relaxing to, need not be scale- invariant [39]. We then 
arrive at the so-called ageing algebra agt(d) := (A 0jl , Y±X 2 , M , Rq )j,k=i,.. d c sd)(d). We 
shall study the consequences of a logarithmic extension of ageing invariance. 

In section 2, we shall write down the generators of logarithmic ageing invariance and shall 
find the co-variant two-point functions in section 3. In section 4, we discuss some applications. 
In particular, we shall show that the scaling of the two-time autoresponse function in ID critical 
directed percolation is well described in terms of logarithmic ageing invariance. We conclude 
in section 5. 



2 Logarithmic extension of the ageing algebra aQt(d) 

For definiteness, we consider the ageing algebra a$je(cf) C sd)(d) of the Schrodinger algebra. 
The generators of the representation (jl.6p can in general be taken over, but with the important 
exception 

77 -4- 1 A/f 77-I-I 
X n = -t n+1 d t — t n r ■ V r - —{n + l)nt n ~ l r 2 — xt n - (77 + l)<t n (2.1) 

where now 77 > and (11. 8p remains valid. In contrast to the representation (II. 6p . we now 
have two distinct scaling dimensions x and £, with important consequences on the form of the 
co- variant two-point functions [HU [36], see also below |§ To simplify the discussion, we shall 
concentrate from now on the temporal part (^f(ti,r)^f(t2,v)), the form of which is described 
by the two generators Xq^, with the commutator [Xi,Xq] = X\. At the end, the spatial part 
is easily added. 

We construct the logarithmic extension of age(rf), analogously to section 1, by considering 
two scaling operators, with both scaling dimensions x and £ identical, and replacing 

in eq. (12. ip . the other generators (II. 6p being kept unchanged. Without restriction of generality, 
one can always achieve either a diagonal form (with x' — 0) or a Jordan form (with x' — 1) of 

4 There is a so-called 'exotic' central extension of CGA(2) 09], but the extension of the known two-point 
functions [3j |H [50] to the logarithmic version has not yet been attempted. 

If one assumes time-translation-invariance, the commutator [X\, X-i] = 2Xo leads to £ = and one is back 

to (rra . 
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the first matrix, but for the moment it is not yet clear if the second matrix in (12. 2p will have 
any particular structure. Setting r = 0, we have from (12.11) the two generators 



Xn = -tdt 



X X 

x 



X 1 = -t 2 d t - t 



x + £ x' + £' 



(2.3) 



and we find [X X ,X Q \ = X 1 + \tx'£" 
must distinguish two cases. 



-1 
1 



= Xi. The condition x'£" = follows and we 



possibilities: (i) either, the matrix 



is diagonalisable. We then 



1. x' = 0. The first matrix in (12 .2\i is diagonal. In this situation, there are two distinct 

e t ) \ o c- 

have a pair of quasi-primary operators, with scaling dimensions (x,£ + ) and (x, £_). This 
reduces to the standard form of non-logarithmic ageing invariance [36]. Or else, (ii), the 

C C \ ^ ( Z i 

situation considered below. 



matrix 



reduces to a Jordan form. This is a special case of the 



2. £" = 0. Both matrices in (12. 2 ft reduce simultaneously to a Jordan form. While one can 
always normalise such that either x' — 1 or else x' — 0, there is no obvious normalisation 
for This is the main case which we shall study in the remainder of this paper. 

In conclusion: without restriction on the generality, we can set £" — in eqs. $2.2Xi2.3\) . 

For illustration and completeness, we give an example of a logarithmically invariant Schrodinger 
equation. Consider the Schrodinger operator 



S := (^Mdt - v; 



2M 



, d 



1 




(2.4) 



Using (I2.3P with the spatial parts restored, we have [S, X ] = —S and [S, X\] = —2tS and 
furthermore, S commutes with all other generators of agz(d). Therefore, the elements of age(c?) 



map any solution of S 



to another solution of the same equation. 



3 Two-point functions 

Consider the following two-point functions, built from the components of quasi-primary oper- 
ators of logarithmic ageing symmetry 



F = F(t h t 2 ) 

G\2 — Gl2\tli ^2) 
G21 = G*2l(tl, £2) 

H = H(t u t 2 ) 



(3.1; 



Their co- variance under the representation (I2.3p . with £" = 0, is expressed by the conditions 
X 0i iF = 0,. . . , where X 0i i stands for the extension of (12.31) to two-body operators. This leads 



5 



to the following system of eight equations for a set of four functions in two variables. 



t ± di + t 2 d 2 + ~{xi + x 2 ) 



t 2 A + t\d 2 + {x x + + (x 2 + 6)^2 F(t 



tidi +t 2 d 2 + - fa +x 2 ) 



G 12 (ti,h) + ^F(t 



t 2 A + t\d 2 + ( Xl + + (x 2 + £ 2 )t 2 G 12 (t u t 2 ) + (x' 2 + Qt 2 F(t 



tidi +t 2 d 2 + ~(xt +x 2 ) 



G 2 i(ti,t 2 ) + ^F(* 1 



tjdi + t\d 2 + ( Xl + + (x 2 + £ 2 )t 2 G 2 i(ti, t 2 ) + K + QhF(t 



t x d x +t 2 d 2 + ~(xi +x 2 ) 



H(t 1 ,t 2 ) + ^G 12 (t u t 2 ) + 3 ^G 21 (t 1 



t 2 ^ + t\d 2 + (x a + + (x 2 + &)t 2 J 
+(ari + ^)tiGi2(ti, t 2 ) + (x' 2 + Qt 2 G 21 (tt 



t2) 











t 2 ) = o 

* 2 ) = (3.2) 

to) = 



t2) 







where di = d/dti. We expect an unique solution, up to normalisations. It is convenient to solve 
the system (13.21) via the ansatz, with y := t\/t 2 



F(t u t 2 ) = t 2 {Xl+X2)/2 yb+te-^V^y - l)-(^i+^)/2-a-6 f( y ) 
G 12 {h,t 2 ) = fr^+^V 2 „&+(*»-*i)/2 



G 21 (t l7 t 2 ) = t -^ +X * )/2 y^-^(y - 1)-^)/^-^^^ . ^.(^ (3.3) 

H(t 1 ,t 2 ) = ^ Xl+X2)/2 



1. The function F does not contain any logarithmic contributions and its scaling function 
satisfies the equation f'(y) = 0, hence 

f(y) = f = cste. (3.4) 

This reproduces the well-known form of non-logarithmic local scaling [36] . 

Comparing this with the usual form (11 .ip of standard LSI with z — 2, the ageing exponents 
a, a', Xr are related to the scaling dimensions as follows: 

1 



a = -(xx+x 2 )-l ' a '~ a = £i + £2 , \r = 2(x 1 + £ 1 ) 



(3.5) 



For example, the exactly solvable ID kinetic Ising model with Glauber dynamics at zero tem- 
perature [25] satisfies (II. ip with the values a = 0, a' — a = — |, Xr — 1, z — 2. Further examples 
of systems with a' — a ^ are given by the non-equilibrium critical dynamics of the kinetic 
Ising model with Glauber dynamics, both for d = 2 and d = 3 
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2. Next, we turn to the function G\ 2 . Co- variance under X n leads to the condition 



(012,1(2/) + \x' 2 f{y) ) + £}0" + 1) in j * 2 • gi2,j+i(y) = o 



which must hold true for all times t 2 . This implies 

1 



0i2,i (y) = -7; x 2f(y) > 9i2,j(y) = ; Vj ^ 0, 1 



In order to simplify the notation for later use, we set 

012(2/) '= 9i2,o (y) , 712(2/) := 012,1(2/) 



-Afiy) 



(3.6) 



(3.7) 



(3.1 



and these two give the only non- vanishing contributions in the ansatz (13. 2j) . Furthermore, the 
last remaining function g 12 is found from the co- variance under Xi, which gives 

ln J h (y{y - l)g' 12jj (y) + (j + l)g w {y)) + (4 + £2)/ (2/) = (3.9) 



3& 



for all times t 2 . Combining the resulting two equations for gi 2 and 712 with ( 13.81) leads to 



y{y - l )g'i 2 {y) 



(3.10) 



3. The function G 2 i is treated similarly. We find 

92i(y) ■■= 92i,o(y) , 721(2/) := 021,1 (y) = 

and the differential equation 



x'J(y) , g 2 iAy) = Q ; for all 0,1 (3.11) 



1 , 



y(y - l) 521 (y) + K + Ci) - -x'J(y) = 



(3.12) 



4. Finally, dilatation-covariance of the function H leads to hj(y) = for all j 7^ 0, 1, 2 and 

1 



h(y) 
h 2 {y) 



{x'i9i2{y) +x' 2 g 21 (y)) 



-x[x' 2 f(y) 



The last remaining function ho(y) is found from co- variance under X\ which leads to 



y(y - l)h' (y) + l(x[ +Qy- -x[j g 12 (y) + i-x' 2 + g 21 (y) = 
Using (|3.4|) . the equations (13. 1Q|I3 . 1 2|I3 . 14j) are readily solved and we find 

012(2/) - ' ' 1 ' " ' * ■' 
021(2/) 
h (y) 



(3.13) 



(3.14) 



012,0+ ( y + £2 ) /oln 



y-i 



021,0 ' 
/i - 



i + ^)/oln|y-l|-^/ o ln|0 



r + f 1 ) 021,0 + ( + ^ ) 012,0 



In \y — 1| — 



x, 



"021,0 



+ £ 2 012,0 



+5/. 



fl + fMlnliz-ll + ^lnl!, 



f +e 2 ) in 2 



ln|j/| 
(3.15) 



7 



where / , <?i 2 ,o, 921,0, h are normalisation constants. We summarise our results: 
F(t 1 ,t 2 ) = t - {xi+X2)/2 y^ + ^- x ^ 2 (y- i)-(*i+*3)/2-*i-fcy 

G 12 (i 1; t 2 ) = ^+*»)/ a j^+^-xV^y - l)-(-X +a;2 )/2-6-6 + ln t2 . 7la(y) ^ 

#(t 1; t 2 ) = t j(*l+ !e »)/ 2 3/ &+( a!a - a »l)/2( ?/ _ 1 )-(xi+«2)/2-{ 1 -& (316) 

x (h (y) + lnt 2 • + In 2 1 2 • ^2(2/) 



where the scaling functions, depending only on y — t 1 /t 2 , are given by eqs. fl3.8f3.llf3.13f3.15p . 

Although the algebra age(d) was written down for a dynamic exponent z — 2, the space- 
independent part of the two-point functions is essentially independent of this feature. The 
change (x,x',£,£ r ) i-> ((2/z)x, (2/z)x', (2/z)£, (2/z)£') in eq. ( 13. 16ft produces the form valid for 
an arbitrary dynamical exponent z. 

Since for z = 2, the space-dependent part of the generators is not affected by the passage 
to the logarithmic theory via the substitution ( 12. 2ft . we recover the same space-dependence as 
for the non-logarithmic theory with z = 2. For example, 



F(ti,tz;r 1 ,r 2 ) = S(M 1 + M 2 ) 9(*x - h) t^ +X2)/2 f 



Mi (n - r 2 ) 

"2 ti - to 



(3.17) 



where we also included the causality condition ti > t 2 , expressed by the Heaviside function 0, 
which can be derived using the methods of [35]. Similar forms hold true for G\ 2 , G21, H. 

Comparison with the result (11.101) of logarithmic Schrodinger-invariance shows: 

1. logarithmic contributions, either as corrections to the scaling behaviour via additional 
powers of ln t 2 , or else in the scaling functions themselves, may be described independently 
in terms of the parameter sets (x' x ,xf^) and (£i,£ 2 )- 

In particular, one may choose to introduce the logarithmic structure only through a single 
one of the two generators X and X\. 

2. If one sets x[ = x' 2 = 0, the scaling functions gx 2 , g 2 \ and h contain logarithmic terms, 
although there is no predicted logarithmic breaking of scaling, in contrast to what occurs 
in logarithmic conformal invariance or logarithmic Schrodinger invariance. 

3. The constraint F = of both logarithmic conformal invariance and logarithmic Schrodinger 
invariance is no longer required. 

4. If time-translation-invariance is assumed, one has £1 = £ 2 = = £ 2 = 0, x\ = x 2 and 
/o = 0. The functional form of eqs. ( I3.16f3.17p then reduces to the Schrodinger-invariant 
forms of eq. (ll.lOp . 
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4 Applications 



4.1 Directed percolation 

It is well-understood that critical 2D percolation can be described in terms of conformal in- 
variance [46]. Notably, Cardy [12] and Watts [72] used conformal invariance to derive their 
celebrate formulae for the crossing probabilities. More recently, it has been shown that a pre- 
cise formulation of the conformal invariance methods required in their derivations actually leads 
to a logarithmic conformal field theory [51] . Since directed percolation is in many respects quite 
analogous to ordinary percolation, we raise the question: 

can one describe dynamical scaling properties of critical directed percolation in terms of 
logarithmic ageing invariance ? 

The directed percolation universality class can be realised in many different ways, with 
often-used examples being either the contact process or else Reggeon field theory, and very 
precise estimates of the location of the critical point and the critical exponents are known, see 
[4TH I5TJ1 155] and references therein, and in agreement with extensive recent experiments [70J. 
In the contact process, a response function can be defined by considering the response of the 
time- dependent particule concentration with respect to a time-dependent particule-production 
rate. The relaxation from an initial state is in many respects quite analogous to what is seen 
in systems with an equilibrium stationary state [T5J EFJ I7J- In figure [TJ we show data of the 
autoresponse function R(t,s) = s -1_a /fl(t/s) of ID critical directed percolation, realised here 
by the contact process. The initial state contains uncorrelated particles at a finite density. The 
data are obtained from the transfer matrix renormalisation group (tmrg) which are consid- 
erably more precise than data obtained from a Monte Carlo simulation, subject to stochastic 
uncertainties [TBI IIH] ■ Aspects of local scaling can be emphasised by plotting the function 

h R (y):=f R (y)y^ z (l-y^Y +a (4.1) 

over against y = t/s, with the exponents taken from [55J. We observe an excellent collapse 
of the data when y is large enough, but we also see that finite-time corrections to dynamical 
scaling arise when y — > 1, the precise form of which depends on the waiting time s. 

Starting from large values of y, and proceeding towards y — > 1, a description of the data in 
terms of local scale-invariance increasingly needs to take finer points into account. First, in its 
most simple form, one would naively assume a = a', when (11. ip predicts a horizontal line in 
this plot. Indeed, this describes the data down to y = t/s « 3 — 4, but fails when y becomes 
smaller. We had tried earlier [36J to take these deviations into account by admitting that a and 
a' can be different. This is equivalent to the assumption that £ + £ ^ and describes the data 
well down to t/s ~ 1.1. However, further systematic deviations exist when t/s is yet closer to 
unity. For the values of s used in figure [TJ it is clear that one is still in the dynamical scaling 
regime and an explanation in terms of a more general form of the scaling function should be 
sought. 

We now try to explain the tmrg data in figure [JJ in terms of logarithmic ageing invariance 
(extended to an arbitrary dynamical exponent z as outlined above). We make the working 
hypothesis R(t, s) = (ip(t)ip(s)), where the two scaling operators ip an d are described by the 
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t/s-1 

Figure 1: Scaling of the autoresponse R(t, s) = s~ 1 ~ a fn(t/s) of the ID critical contact process, 
as a function of y = t/s, for several values of the waiting time s, and indicated by the dash- 
dotted lines. The dashed line labelled 'lsi' gives the prediction for the scaling function hn(y) = 
fn(y)y XR l z (l — l/y) 1+a as obtained from standard, non logarithmic local scale-invariance (11.11) . 
with a' — a = 0.26. The full curve labelled 'LSI loga' is the prediction (14. 3f) of logarithmic local 
scale-invariance with £' = 0, as described in the text. 



logarithmically extended scaling dimensions 

In principle, one might have logarithmic corrections to scaling, according to eq. (I3.16p . Because 
of the excellent scaling behaviour seen in figure (TJ we conclude that logarithmic corrections are 
absent in the data. Hence the two functions hi^{y) must vanish. Because of eq. (I3.13p . this 
means that x' = x' = 0. Then logarithmic ageing invariance (I3.15P predicts 

h R (y) = (l-^j ^o-torin(l-l/y)-^o? 2 ln 2 (l-l/y) 

- 92id'Hy - 1) + \M' 2 ln 2 (2/ - 1)) (4.3) 

Since for y sufficiently large, the numerically observed scaling function h^{y) becomes essen- 
tially constant, we conclude that there are no leading logarithmic contributions in the y — > oo 
asymptotic behaviour, hence £' = and the second line in f !4.3p vanishes. Hence we arrive at the 
following phenomenological form hn(y) = h (l — 1 / y) a ~ a ' (l — v41n(l — 1/y) — I?ln 2 (l — l/y)), 
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with the normalisation constant ho and where the universal parameters A, B and the exponent 
a — a' must be determined from the data. The full curve in figure [I] shows to what extent this 
describes the available data, with the chosen values 



a' - a ~ 0.174 , A ~ 0.13 , B ~ 



0.0168 



— 



0.0888 



(4.4) 



Indeed, the chosen form gives a full account of the tmrg data, down to t/s — 1 ~ 2 • 10~ 3 , 
about two orders of magnitude smaller than for non-logarithmic local scale-invariance. This is 
about the size of the region where dynamical scaling has been confirmed through the collapse 
of data with different values of s. However, we also point out that the functional form of hn{y) 
may depend quite sensitively on the values of the several parameters such that error bars in 
a' — a, A,B are of the order of at least 25%. 

In particular, the tmrg data suggest that the logarithmic nature of the scaling operator 
should merely enter via the response field, and only through the consideration of the 'special' 
transformation X%, since x' = x = £ = and £' 7^ is the only quantity which describes the 
departure from the standard non-logarithmic scaling. We observe that the present estimate 
a' — a = 0.17(5) is considerably more small than our earlier estimate a' — a w 0.27. This is 
the first time that a theory could be formulated which describes the autoresponse in the entire 
range of the scaling variable, 2 • 10~ 3 < y — 1 < 00. We point out that existing field-theoretical 
methods based on the ^-expansion [7J obtain reliable results only in the opposite case y ^> 1, 
notably on non-equilibrium exponents and universal amplitudes. 

4.2 Logarithmic scaling forms 

In the ageing of several magnetic systems, such as the 2D XY model quenched from a fully 
disordered initial state to a temperature T < Tkt below the Kosterlitz-Thouless transition 
temperature [HI |9l [1] or fully frustrated spin systems quenched onto their critical point [TTJ HI] , 
the following phenomenological scaling behaviour 



has been found to describe the simulational data well. Could this scaling form be explained 
within the context of logarithmic ageing invariance ? Helas, this question has to be answered 
in the negative. If one fixes y = t/s and expands the quotient lns/lnt = lns/(ln?/ + Ins) for 
s — > 00, eq. (14. 5 p leads to the following generic scaling behaviour 



Comparsion with the explicit scaling forms derived in section 3 shows that there arise only 
combinations of the form ln n y ■ ln m s or \n n (y — 1) • In" 1 s, where the integers n, m must satisfy 
< n + m < 2. This is incompatible with (I4.6p . 

In conclusion, the logarithmic scaling form (14. 5p cannot be understood in terms of logarith- 
mic ageing invariance, as presently formulated. 




(4.5) 




(4.6) 
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5 Conclusions 



We have discussed the extension of dynamical scaling towards local scale-invariance in the case 
when the physical scaling operator acquires a partner with the same scaling dimension. Since in 
far-from-equilibrium relaxation, time-translation-invariance does not hold, one cannot appeal 
directly to the known cases of logarithmic conformal and Schrodinger-invariance. Indeed, anal- 
ogously to the non-logarithmic case, the doubletts of scaling oprators are described by pairs of 
Jordan matrices of scaling dimensions. When computing the co- variant two-point functions, the 
absence of time-translation-invariance allows, independently, to include logarithmic corrections 
to scaling and also non-trivial modification of the scaling functions, see eq. fl3.15f3.16p . This 
generalises the forms found from logarithmic conformal or Schrodinger-invariance [JT]. 

Motivated by the fact that imporant properties of ordinary 2D critical percolation can be 
understood in terms of logarithmic conformal invariance [5T] . we have re-analysed the autore- 
sponse R(t, s) of critical ID directed percolation in terms of our logarithmic extension of local 
scale-invariance. The available data suggest that at least this observable behaves as if directed 
percolation were described by logarithmic local scale-invariance (with an obvious generalisation 
to z 7^ 2). Of course, further independent tests of this possibility are required. 

Since logarithmic conformal invariance also arises in disordered systems at equilibrium, it 
would be of interest to see whether logarithmic local scale-invariance could help in improving 
the understanding of the relaxation processes of disordered systems far from equilibrium, see 
e.g. [6311371 SHIES]. 

Acknowledgement: I thank T. Enss for the tmrg data, M. Pleimling for useful correspon- 
dence and the Departamento de Fisica da Universidade de Aveiro for warm hospitality. 
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